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ABSTRACT 


The  terminal  ballistic  dispersion  of  a  non-rotational,  small-caliber 
weapon  system  is  obtained  on  a  priori  knowledge  of  the  initial  dispersion 
induced  by  the  weapon  syicem  proper. 

The  study  considers  a  dynamic  system  whose  behavior  is  governed  by 
a  system  of  differencial  equations  having  probabilistic  initial  conditions. 
The  behavior  of  the  system,  in  terms  of  the  kinematic  probability,  is  then 
determined  as  a  function  of  time  and  spatial  variables. 

Subsequently,  in  contrast  to  the  customary  method  in  the  evaluation  of 
weapon  systea  effectiveness,  a  measure  of  effectiveness  -  probability  of 
hit  -  is  obtained  as  a  function  of  initial  dispersion.  Also,  numerical 
examples  as  well  as  discussion  of  the  results  are  given. 


I.  INTRODUCTION 


Conventionally,  the  terminal  effect  of  projectiles,  related  to  a 
weapon  system,  is  analyzed  based  on  a  priori  knowledge  of  the  dispersion 
at  a  target.  The  probability  density  function  of  the  dispersion  is 
specified  and  its  parameters  are  given.  The  probability  of  hit  and  other 
effectiveness  measures  are  then  sought  [1,2]. 

Such  an  approach  is  pertinent  if  we  address  ourselves  strictly  to  the 
problem  of  effectiveness  at  the  target  only,  and  the  terminal  dispersion 
taken  is  well  substantiated  experimentally.  When  the  projectiles  are 
related  to  a  new  conceptual  weapon  system  for  which  no  physical  experiments 
have  yet  been  performed,  a  question  naturally  arises  about  the  soundness 
of  the  above  approach. 

To  a  larger  extent  than  this  consideration  on  the  target  alone,  if 
we  are  more  concerned  with  the  effectiveness  of  a  conceptual  weapon  system 
which  incurs  certain  terminal  effects  on  a  given  target,  hoping  eventually 
to  bolster  the  rationale  for  preliminary  engineering  design  of  the  system, 
it  is  hardly  plausible  that  we  could  merely  consider  the  terminal  dispersion 
as  a  priori  knowledge.  What  appears  to  be  needed  in  the  treatment  of  such 
a  problem  is  a  terminal  dispersion  obtained  on  the  base  of  some  criteria 
pertaining  to  the  weapon  system  proper. 

We  can  schematically  envision  three  regions  of  concern:  conceptual 
weapon  system,  exterior  ballistics  and  terminal  effects,  as  shown  in  Fig.  1. 
The  region  D  which  is  the  intersection  of  the  regions  A  and  B  contains 
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information  relating  to  the  weapon  system,  whereas  the  region  E  whiv.h  is 
the  intersection  01  the  regions  B  and  C  has  information  pertaining  to  the 
terminal  dispersion.  It  is  clear  that  information  arises  from  D,  coupled 
with  B,  would  arrect  information  about  E  and  subsequently  about  C.  The 
specific  point  of  concern  now  is  what  particular  information  from  D  must 
we  generate  and  how  should  it  arrive  at  E.' 

Since  a  measurable  Information  at  E  is  the  terminal  dispersion,  it 
Is  reasonable  to  consider  its  counterpart  at  D  -  the  initial  dispersion 
For  this  reason,  the  a  prion  knowledge  on  the  Initial  dispersion  of  a 
conceptual  weapon  system  is  necessary,  for  it  is  formulated  on  the  base  o' 
our  understanding  of  the  intrinsic  properties  such  as  weapon  dynamic 
parameters  [3]  ot  the  system  under  consideration  The  terminal  dispersion 
is  then  rigorously  sought  analytically 

In  this  report  the  probability  distribution  or  terminal  dispersion 
and  subsequently  some  terminal  eirects  are  obtained  lot  a  given  target  and 
a  given  probability  ot  initial  dispersion  Tne  problem  is  considered  as  a 
stochastic  dynamic  process  tor  a  time  duration  ranging  from  t,  u  t 
The  initial  probabilities  of  the  kinematic,  or  state,  variables  -  displace¬ 
ment  and  velocity  -  are  given  a  prion- 

As  of  special  interest,  an  exterior  ballistic  model  tor  small  mass  is 
considered.  It  is  shown  that  a  linearization  or  the  system  dynamic  equation 
can  be  realized  A  method  of  attack  rendering  the  transformation  ot  the 
kinematic  probabilities  from  one  time  to  another  is  described  Solutions 
of  the  problem  are  then  presented  A  numerical  example  is  given,  and  the 
results  are  discussed. 
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II.  STATEMENT  OF  THE  PROBLEM 

Let  us  consider  a  weapon  ayscem  of  caliber  d  delivering  a  small 
projectile  of  mass  m  into  a  target.  The  projectile  has  a  muzzle  velo¬ 
city,  V  ,  with  respect  to  a  given  rectangular  Cartesian  coordinate  system, 
x-y-z,  which  is  also  an  inertia  frame,  having  the  y  coordinate  axis  in 
the  range  direction.  The  range  to  the  target  is  short  such  that  the  air 
density  p  and  the  temperature  in  the  neighborhood  of  the  trajectory,  as 
well  as  the  gravitational  acceleration,  g,  can  be  taken  as  constant.  It 
follows  the  velocity  of  sound  aQ  is  also  constant. 

Also  let  us  assume  the  air  is  almost  still  and  the  projectile  moves 
with  its  axis  tangent  to  its  trajectory.  The  only  force  acting  on  the 
projectile  are  the  drag  and  the  gravitational  force.  We  can  now  write 
down  the  normal  equations  of  motion  for  the  projectile  as  follows  [4] 

x  -  -KjjPd2  V  V^/m 

y  -  -l^pd2  V  Vy/m  (2.1) 

z  -  -KDpd2  V  V^/m  -  g 

where  V  ,V  and  V  are  the  components  of  V  along  x,y  and  z  coordinate 
x  y  2 

axes  respectively,  and  is  the  drag  coefficient. 

It  has  been  shown  in  [5,6,7]  that  can  be  considered  as  a  function 

of  Mach  number  for  (2.1)  and  for  certain  range  of  Mach  numbers  is 

inversely  proportional  to  the  Mach  number  of  the  projectile.  Consequently, 
it  can  be  written  as  follows 
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where  6  is  a  positive  constant  cf  proportionality  and  depends  on  the 
type  of  projectile,  in  particular,  the  geometrical  configuration. 

A  simple  linearization  of  (2.1)  can  be  done  by  substituting  (2.2) 
into  (2.1). 

x  ■  yx 

y  -  yy  (2.3) 

z  •  yz  -  g 

where  y  -  -0ao  p  d2/m. 

Letting  x  -  x,x2  -  y  and  x^  -  z,  we  can  reduce  (2.3)  to  a  system 
of  first  order  linear  system,  with  the  initial  conditions  at  tQ  ■  0 
specified,  as  follows: 

*i  -  xi+3  i-1,2,3 

x^  «  yx^  j-4,5  (2.4) 

*6  "  YVg 

and  xi(0)  ■  x^  i*l,2,  •••,6 

Now,  x°  is  not  deterministic  for  there  exists  initial  dispersion. 
Therefore,  (2.4)  is  a  system  of  stochastic  differential  equations  with 
stochastic  initial  conditions. 

Our  problem  can  now  be  stated  succinctly  as  follows:  given  the  initial 
Joint  probability  density  function  of  x°,  find  the  joint  and  marginal  pro¬ 
bability  density  functions  of  x^(t)  for  (2.4). 
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A  natural  follow-up  question  is  as  follows:  given  a  target  T,  find 
the  probability  of  x^(t)  over  T,  which  is  the  probability  of  hit. 

In  this  report,  we  presume  that  x°,  a  random  vector  with  components 
x° ,  has  normal  joint  probability  density  function  with  independent  marginal 
probability  density  functions  characterized  by  two  parameters  -  mean  i.? 
and  variance  (c°)2.  A  rectangular  parallelepiped  target,  as  shown  in 
Figure  2,  is  considered,  i.e, 

T  -  {(x  ,x  ,x  )eE3|s  <_  x  <_  s  ,  s  <  x  <  s  ,  s  <  x  <  s  >  (2.5) 

l  i  —  l  3  2  ^  5  3  6 

where  s^  is  specified. 
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Ill-  METHOD  OF  SOLUTION 


For  convenience  let  us  consider,  in  general,  a  system  of  n  first 
order  duterential  equations  written  in  matrix  form 


y  -  £(y,t)  t  -  0 

y(0)  -  yL 

where  ^ 

y  ■  (y  »y  ,  •••  ,y„) 

i  n 

f  -  (i  ,f,,  •••  ,tn)T 

and 

,00  .T 

y(0)  -  (y  ,y2,  .yn>  • 


(3.1) 


Geometrically  (3.1)  describes  the  dynamics  of  a  point  in  En  with 
time  t  as  a  parameter.  If  the  initial  condition  is  known  only  probabi¬ 
listically,  then  the  solution  of  the  vector  equation  (3.1)  is  a  random 
vector  even  though  the  equation  itself  is  deterministic.  Thus,  we  can 
solve  for  the  solution  y,  a  random  vector,  by  considering  these  equations 
as  if  they  were  deterministic  ones,  and  obtain  the  solution  in  the  form 

y=*(y\c)  (3.2) 


We  observe  that  in  (3.2)  #  is  a  deterministic  transformation  mapping  the 

random  vector  y°  denned  on  a  probability  space  ft  into  a  new  random 
vector  y  derined  on  the  same  space  as  shown  in  Figure  2A.  Moreover,  this 
mapping  is  deterministic.  Hence,  <j>  is  probability  -  preserving  i.e. 

P(y^A)  =  P(>°eA)  (3,3) 


where  P  denotes  a  probability  distribution  and  A  is  a  measurable  cibse: 


of  the  reals- 
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These  notions  lead  us  readily  to  an  approach  for  solving  (2.4)  based 
on  the  previous  works  [8,9].  For  convenience,  let  us  recapitulate  the 
results  in  a  form  of  a  theorem  without  proof: 


Theorem 


Given  the  system  (3.1)  and  the  initial  Joint  probability  density 
function  p(y°)  at  t  -  t  .  If  for  all  i,j-l,2,  •••  ,n 


3fi 

t~~~  I  <  M 

9yj  - 


(3.4) 


in  the  domain  of  definition,  where  M  is  a  positive  constant.  Then,  the 
probability  density  function  for  t  >  tQ,  denoted  by  p(y,t),  which 
satisfies  the  following  partial  differential  equation 


n 

l 

i-1 


3(pf1) 

3yi 


0  , 


(3.5) 


is  given  by  the  formula 

p(y,t)  -  p(y°)  |  |  (3.6) 

where  |^—  is  the  Jacobian  of  the  initial  vector  y°  with  respect  to 
the  vector  y. 

Note  that  it  is  essential  for  this  approach  that  the  inverse  of  y 
in  (3.2),  i.e., 

y0  -  *(y)  (3.7) 

where  .  -1 

*  ■  4>  » 

must  be  attainable  in  order  to  utilize  (3.6). 
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IV.  SOLUTION  OF  THE  PROBLEM 

We  observe  that  the  theorem  is  applicable  to  (2.4).  Hence,  our 
problem  reduces  to  obtain  the  inverse  of  the  solution  of  (2.4). 

Let  us  denote 


/  0  0  0  \T 

(VV  •••  •x6) 


and 


(x  ,x  ,  •• • 

1  2 


|X  )  • 
6 


Solving  (2.4)  for  x  ,  we  obtain 
x3  -  Ax+B 


where 


0 

1 

0 

0 

0 

0 


0 

0 

1 

0 

0 

0 


a 

0 

0 

b 

0 

0 


0 

a 

0 

0 

b 

0 


and 


B  -  (0,  0,  ~*<**t}  ,  0,  0,  -ag)T 


a  =  (b-l)/y 


b  -  e 


-Yt 


The  absolute  value  of  the  Jatobian  of  (4.2)  is  then 

I  lr  i  -  e"3n 


(4.1) 

(4.2) 

(4.4) 

(4.5) 

(4.6) 

(4.7) 
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Now,  Che  joint  probability  density  function  of  x  at  any  time  t  can 
be  written  by  using  (4.2)  with  the  normal  joint  probability  density  function 
of  x°  as  follows: 


-3yt  ,  6  x  -v 

p(x,t)  -  - 7 -  exp  {-  J  l  (  — - —  )  ) 

/0  ,3  0  i-1  oj 

(2 u)  no.  i 

i-1  1 

where 

^  m  uj  i-1, 2, 4, 5 

V3  •  Ul  +  y  (t+a) 
and  Q 

v6  -  u6  +  ag. 


(4.8) 


We  have  tacitly  used  the  independence  among  the  marginal  probability 
densities  of  x°  as  stated. 

It  follows  that  the  marginal  probability  density  of  x  is  expressible 


as 


where 


p(x  ,t)  -  I  ,  p(x,t)dx 
-  E 


dr  -  tt  dx . 

J-l  J 
J*i 


for  n-6,  and 


(4.9) 


the  integral  sign  denotes  n-1  integrals  over  En  1.  Thus  we  can  write 
with  the  use  of  (4.8)  and  (4.9)  following  results 


p(*1»t)  - 


exp{-  y  [  — 


x.  -  <y»  2 

n  ■*  T 


(4.10 


where 


8i  “  (a  °i+3)2  +  (b  °i)2 


i-1, 2, 3 
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and 


p(x.,t)  - 


/27 


f  1  /  bXi"Vi  .2. 
exp  [-  y  (  - 5 —  )  1 


i-4,5,6 


(4.11) 


The  mean  of  is  expressible  by  definition  as 

r  ® 

u1(t)  »  I  xt  p(x1,t)dx1  (4.12) 

which  yields  the  following  results  with  the  use  of  (4.10)  and  (4.11): 


ut(t)  -  ^  -  avi+3/b, 

>1,2,3 

(4.13) 

^(t)  -  v1/b 

i-4,5,6 

(4.14) 

Similarly,  the  mean  square  value  of 

f  °° 

E[x^(t)J  -  |  x2  p(xi,t)dxi 

*  -X 

which  gives 

x^  can  be  written  as 

(4.15) 

0,  av.  - 

ElxJU)]  -  )  -  O,  ~  l 

0 

)2,  i-1.2,3 

(4.16) 

E[x2 (t) ]  -  (oj  f  v2)/b2 

i-4,5,6 

(4.17) 

Since  the  variance  of  x^  is  related  to  its  mean  and  mean  square  value  as 

o2(t)  -  E[x2(t)]  -  u2(t)  (4.18) 

the  following  results  are  immediate 


•  * 

: 
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i-1,2,3 


(4.19) 


3*(t)  -  81/b2 
o2( t)  -  [o°/b]2 


i-4,5,6 


(4.20) 


For  che  target  T  of  (2.5)  the  probability  of  hit  at  a  given  time  t, 

P,  (t),  can  be  obtained  by  firat  finding  the  joint  probability  density 
n 

function  of  x^,  i*l,2,3;  i.e. 

f 

P (x  ,x  ,x  ,t)  -  I  p(x,t)dlf  (4.21) 

'E 

where 

dx  -  dx  dxc  dx, 

«  5  6 


Subsequently,  P^(t)  over  T  can  be  written  as  follows 


ph(t) 


p(xj ,x2,x3,t)dx 


(4.22) 


where  dx  ■  dx  dx2  dx  .  and  the  integral  u.'yr  denotes  a  triple  integral 
over  T. 

Substantially,  by  utilizing  the  results  of  (4.21)  into  (4.22)  and 
employing  the  notion  of  error  functions,  we  obtain  the  following  expression:' 
3 

ph(t)  ■  n  [erf(n2i)  -  erf (n21_1)]/2  (4.23) 

where 


n2i-l  "  b(82i-l  "  lJi)/'/*®i 


n2i  "  b(s2i  “  ui‘) / 
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V.  NUMERICAL  EXAMPLE 


As  an  illustration,  let  us  consider 

the  following  weapon  system  in 

the  milieu  as  given: 

-4 

m  -  3. 774x10 

slug  (85  grain) 

-2 

d  -  1.9685x10 

ft  (6.0mm) 

8  -  0.25 

aQ  -  1120.27 

ft/sec 

-3 

3 

p  -  2.377x10 

slug/ ft 

g  -  32.174 

f t/sec^ 

The  system  gives  initial  dispersion  with  means  (ft)  and  standard 

deviations  (ft)  as  follows: 

u°  -  3.xl0'5 

1 

o°  -  1.5x10"’* 

i 

_°  -  -6xl0"5 

c°  -  3.xl0‘5 

l  °  ■  2x10  ^ 

o  j  -  10'5 

^  u  “  2.0  (ft/ sec) 

o°  »  0.1  (ft/sec) 

■  3500.0  (rt/sec) 

o°  -  50.0  (ft/sec) 

u°,  ■  3.0  (ft/sec) 

b 

o°  -  0.2  (ft/sec) 

A  rectangular  parallelepiped  target 

T  has  the  following  dimensions 

(ft)  for  its  s  : 

s  -  -3.0 

s2  -  3.0 

s,  -  -2.0 

O 

CM 

■ 

• 

s5  -  -3.5 

s,  -  3.5 

b 

12 


f 


Let  ua  consider  Chat  T  la  zaroad  in,  and  la  situated  at  a  range  x" . 
Since  T  la  a  volume  target,  let  us  specifically  address  ourselves  to  a 
particular  time  t1  which  la  related  to  x’  with  nil  initial  variance. 

A  deterministic  relation  can  be  readily  obtained  for  t1  from  (2.4)  as 
follows: 

t1  -  ^  ln[l  +  y(*2  "  *2^*5^ 

As  an  illustration,  a  set  of  ranges  up  to  x‘2  -  1000  meters  is  con¬ 
sidered.  For  any  value  of  x*  in  this  set,  (5.1)  yields  the  corresponding 
t1,  and  the  results  in  Section  4  are  then  subaequently  utilized  in  seelcing 
numerical  solutions. 

In  Fig.  3  the  "flight  time"  t‘  across  the  range  x2  of  interest  is 
given.  The  rengewerd  mean  velocity  w5  is  shown  in  Fig.  4.  In  Fig.  5 
the  mean  lateral  displacement  and  velocity  u  and  are  shown.  The 
rangeward  mean  transverse  displacement  and  velocity,  u-3  and  wQ  are  given 
in  Fig.  6. 

Fig.  7,  8  and  9  depict  the  rangeward  standard  deviations,  and 
°i+3'  r«P«ctively-  Th*  maxima  of  Che  marginal  density  function  p(x^) 
ere  given  in  Fig.  10,  11  and  12.  Finally,  the  effects  of  target  dimensions 
and  initial  dispersion  on  the  probability  of  hit  P^tt  )  are  shown, 
respectively,  in  Fig.  13  and  14. 
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VI.  DISCUSSION  AND  CONCLUSIONS 


A  solucion  to  the  problem  of  finding  the  kinematic  probability  of 
terminal  ballistics,  given  an  initial  dispersion,  has  been  presented  in 
this  report.  The  feasibility  of  such  a  solution  implies  that  the  initial 
dispersion  of  a  weapon  system,  rather  than  the  terminal  one,  should  be 
addressed  rigorously  in  a  weapon  system  effectiveness  study. 

For  any  given  range  and  target  configuration,  the  terminal  effective¬ 
ness,  such  as  probability  of  hit,  is  simply  a  consequential  result  once 
the  initial  dispersion  of  the  weapon  proper  is  attained.  In  contriving 
the  preliminary  design  layout  of  a  conceptual  weapon  system,  the  initial 
dispersion  is  therefore  one  of  the  most  significant  factors  for  considera¬ 
tions.  It  is  a  key  juncture  to  the  overall  system  effectiveness. 

The  method  of  attack  presented  here  is  general  to  the  extent  that  no 
constraint  is  posed  on  the  probability  of  initial  dispersion.  Indeed,  the 
methods  can  also  be  employed  for  the  cases  other  than  normal  distribution 
as  specified  for  the  problem.  The  theorem  is  applicable  to  a  nonlinear 
system  as  well  as  a  linear  one  whenever  (3.7)  is  obtainable. 

As  we  can  observe  from  (4.2)  and  (4.3),  hence  X°  are  coupled 

among  each  other.  The  common  concern  for  an  Initial  dispersion  in  terms 
of  the  initial  displacement  alone  is  therefore  not  adequate.  Clearly,  we 
should  also  take  the  initial  velocity  due  to  weapon  dynamics  in  addition  to 
the  muzzle  velocity,  into  account  for  consideration  of  an  initial  dispersion. 

The  result  of  (4,23)  for  a  volume  target  can  readily  be  applied  to  an 
area  target  of  rectangular  shape  by  collapsing  one  dimension  and  letting 
one-half  of  the  difference  of  the  error  functions  for  that  coordinate 


assume  cha  value  of  unity  In  computation.  Moreover,  (4.23)  indicates 

that  P.  (t)  for  a  volume  target  is  bounded  by  one-half  of  the  corresponding 
n 

value  for  the  area  targec. 

Using  (4.14)  the  absolute  value  of  the  mean  striking  velocity  ^v(t) 
can  be  written  as 

6  1/2 

u  (t)  -  [  Z  V  r'Vb  (6.1) 

i-4 


We  observe  as  uv(t)  vanishes  as  t  increases  to  infinity.  On  the  other 


hand,  from  (4.13)  we  can  write  the  asymptotic  value  for  the  mean  dis¬ 
placement  ui(t),  1-1,2 ,3  as  follows 


Vl<0  ■*  + 


i+3 

TTf 


(6.2) 


which  provides  an  upper  bound  tor  terminal  mean  displacement.  The  para¬ 
meter  y  In  (2.4)  and  (4.6)  can  lead  us  to  consider  1/Iy|  ns  the  time 
constant  of  the  system. 

The  direction  cosines  of  the  striking  velocity  vector  are 
obtainable  from  the  ratio  of  (4.14)  and  (6.1)  as  follows: 


cos  a.  -  — z 
x  6 


i-4, 5, 6 


(6.3) 


I  I  v2)1'2 

j-4  2 


We  observe  that  is  time  invariant;  i.e.,  the  at  the  terminal  point 

coincides  with  that  of  any  departure  point. 

It  should  be  noted  that  the  modeling  does  not  Involve  rotational 
effect  of  a  projectile  Furthermore,  the  relation  (2.2)  should  be  observed. 
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The  advantage  in  choosing  in  the  fora  (2.2)  is  two-fold.  It  not 

only  reduces  the  nonlinear  systems  of  (2.1)  into  a  linear  system  but 
also  gives  a  "better"  approximation  to  the  empirical  curves  of  for 

the  ranges  of  our  concern  (around  1  to  5  Machs),  in  contrast  to  the 
customary  assumption  being  constant. 

The  numerical  results  in  Fig.  4  through  6  can  be  obtained  identically 
by  solving  (2.4)  as  a  deterministic  system  with  x°  -  The  standard 

deviations  shown  in  Fig.  7  to  9  are  almost  linear  with  respect  to  the 
range.  The  maximum  of  displacement  probability  density  function, 
p(x^ )  ,  i-1, 2, 3,  decays  rapidly  before  400m  as  shown  in  Fig.  10,  11  and 

12.  In  the  same  figures,  the  maximum  of  velocity  probability  density 
function  is  shown  having  slow  diversion  Initially  up  to  a  point  about  half 
way  of  the  range,  and  then  an  escalation. 

In  Fig.  13,  it  la  shown  that  as  the  target  dimension  Increases  the 
probability  of  hit  P.^Ct  1 )  Increases  as  well.  When  the  initial  dispersion 
decreases,  P^(t‘)  increases  as  shown  in  Fig.  14.  The  sensitivities  of 
p(x^)  and  P^(t‘)  with  respect  to  system  parameters  and  initial  dispersion 
can  be  further  investigated  by  using  the  results  of  the  Section  4. 
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Effectiveness  Sttidv  of  a  Conceptual  Weapon  System 


Fig.  2  E  Space  With  Coordinates  and 


Function  Space  of  all  Random 
Variables  on  ft 
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Fig.  2A  Schematic  Representation  of  a  Measure 
Preserving  Transformation  d> 


2.0 


* 


u 

<D 

« 


1.5  H 


1.0  t- 


0.5  U 


0.0 


0  200  400  600  800  1000 

Range  x2  meters 

Fig.  5  Mean  Lateral  Displacement  &  Velocity  V9.  Range 


21 


FI*.  14  Effect  of  Initial  Dispersion  on  Probability  of  Hit 
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llity  Density  Functions  vs.  Range 
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Fig.  12  Maximum  of  Probability  Density  Functions  vs.  Range 


Fig.  13  Effects  of  Target  Dimensions  on  Probability  of  Hit 
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